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Bose-Einstein correlations of two identically charged Q- 
bosons are derived considering these particles to be confined 
in finite volumes. Boundary effects on single Q-boson spec¬ 
trum are also studied. We illustrate the effects on the spec¬ 
trum and on the two-Q-boson correlation function by means 
of two toy models. We also derive a generalized expression for 
the Wigner function depending on the deformation parameter 
Q, which is reduced to its original functional form in the limit 
of Q —> 1. 


I. INTRODUCTION 


Two-boson interferometry has for a long time been 
linked to high energy heavy-ion collisions as one of the 
tools to probe the existence of a new phase of matter 
of strongly interaction particles, the quark-gluon plasma 
(QGP), at high temperature and high baryon density 
[1,2]. The hope of discovering the QGP in high en¬ 
ergy heavy-ion collisions is to some extent connected to 
the possibility of measuring the geometrical sizes of the 
emission region of secondary particles. And that is the 
connection point, i.e., so-called Hanbury-Brown-Twiss 
(HBT) interferometry [3,4] method, originally proposed 
in the 50’s for measuring stellar radii. This method has 
been largely studied over the last twenty years, and has 
extensively been developed and improved ever since [5]. 

In a previous paper [6], we have studied boundary ef¬ 
fects on the single-particle distribution and on the two- 
particle correlation function, motivated by the need to 
consider more realistic finite systems, and by the idea 
suggested in Ref. [7]. In that reference, it was shown that 
in heavy-ion collisions the pion system could be thought 
as a liquid of quasi-pions subjected to a surface tension. 
Naturally, it would be expected that this surface tension 
would affect the spectrum distribution, which was shown 
in Ref. [6,8-12]. As pion interferometry is sensitive to the 
geometrical size of the emission region as well as to the 
underlying dynamics, we would expect that the bound¬ 
ary would also affect the correlation function, which was 
indeed demonstrated in Ref. [6]. 

Some time ago, on the other hand, the concept of quons 
was suggested [13] in association to a deformation pa¬ 
rameter, Q, which was viewed as an effective parameter 


able to encapsulate many essential features of complex 
dynamics of different systems. (We call the attention to 
the notation adopted throughout the paper. We use cap¬ 
ital Q to refer to the bosons under study here for avoid¬ 
ing confusion with the relative momentum of the bosonic 
pairs, q = Pi — p 2 , commonly used in interferometry and 
adopted here as well). Effectively, the way it works is 
by reducing the complexity of the interacting systems 
under study into simpler relations, nevertheless at the 
expense of deforming their commutation relations, and 
thus making these more complicated. This is known as 
Q-deformed algebras, an approach which has been widely 
studied in statistical physics [14] and also in heavy-ion 
collisions [15]. Particularly interesting is the approach 
in Ref. [16], where it was shown that the composite na¬ 
ture of the particles (pseudo-scalar mesons) under study 
could result into Q-deformed structures linked to the de¬ 
formation parameter Q. In that reference this parameter 
is then interpreted as a measure of effects coming from 
the internal degrees of freedom of composite particles 
(mesons, in our case), being the value of Q dependent on 
the degree of overlap of the extended structure of the par¬ 
ticles in the medium. Being so, the Q-parameter could 
be related to the power of probing lenses, for mimicking 
the effects of internal constituents of the bosons. In this 
case, and for high enough magnification, the bosonic be¬ 
havior of the Q-bosons could be blurred by the fermionic 
effect of their internal constituents, which would result 
in decreasing the value of Q. We will see that our results 
are also compatible with this interpretation. 

In view of our previous study of confined pions sub¬ 
jected to finite size boundaries, and of the Q-boson ap¬ 
proach mentioned above, we realized it would be very 
interesting to analyze its effects on the two-identical Q- 
boson correlation function. Besides, adding this extra 
degree of freedom extends and generalizes our previous 
approach. Along these lines. Ref. [15,17] turned out to be 
of special interest to our investigation. However, in those 
references the approach was focused on the intercept, A, 
of the two-particle correlation function at zero momen¬ 
tum difference, (i.e., A = C(q — 0, K) — 1), and restricted 
to single modes only. All the possible consequences on 
the effective geometrical information, which are gener¬ 
ally even more interesting, were completely neglected. In 
this paper we develop full Q-boson two-particle interfero¬ 
metric relations and simultaneously study the additional 
effects caused by the finite size boundary on the Q-boson 
spectrum and on the two-Q-boson correlation function. 


I 


The plan of this paper is as follows: in section II, we 
derive the Q"boson single-inclusive distribution, as well 
as the two-Q-boson correlation function, considering a 
density matrix suited for describing charged identical Q- 
boson correlation effects. In section III, the boundary 
effects on the two-Q-boson correlation and single parti¬ 
cle spectrum distribution are illustrated by means of two 
simple specific examples. The conclusions are discussed 
in section IV. Finally, we discuss two complementary top¬ 
ics in Appendices I and II. In the hrst one, we discuss the 
limit Q —> 0 in detail. In Appendix II, we re-derive the 
relations for the single- and two-particle distributions, as 
well as the generalized version of the Wigner function for 
another type of Q—boson [15,17], different than the one 
discussed in the body of the paper. 

II. SPECTRUM AND TWO-Q-BOSON 
CORRELATION FUNCTION 

In this section, we derive general relations for describ¬ 
ing the single- as well as the two-particle inclusive dis¬ 
tributions, which would be suited for describing charged 
Q-bosons bounded in a hnite volume. For doing this, we 
extend the hypotheses assumed in Ref. [6] to the pions 
considered here as a Q-boson system. Essentially, these 
could be summarized as follows: the effects of interac¬ 
tion among the Q-bosons could be modeled by consid¬ 
ering that they move in an attractive mean field poten¬ 
tial, which extends over the whole system. In the two- 
(quasi)particle case, this implies that they would not suf¬ 
fer any other effects besides the mean field attraction and 
the identical particle symmetrization. The effect due to 
the fermionic (constituents) internal degrees of freedom, 
along the lines suggested in Ref. [16], if any, would be 
represented by the effective deformation parameter Q. In 
the present analysis, as assumed before in Ref. [6], the pi¬ 
ons represented by Q-bosons are considered to be quasi¬ 
bound in the system, with the surface tension [7] acting as 
a reflecting boundary. The Q-boson wave function could 
then be considered as vanishing outside this boundary. 
Once more, we assume that these particles become free 
when their average separation is larger than their interac¬ 
tion range and we consider this transition to happen very 
rapidly, in such a way that the momentum distribution 
of the Q-bosons could be essentially governed by their 
momentum distribution just before they freeze out. We 
then study the modifications on the observed Q'boson 
momentum distribution caused by the presence of this 
boundary. On the other hand, we know that interferom¬ 
etry is sensitive to the geometrical size of the emission 
region as well as to the underlying dynamics, and, we 
would expect that the boundary would also affect the 
correlation function, similarly as it affected the pions in 
[6]. However, as we shall see later, there is a significant 
difference in the present case: the parameter A, i.e., the 


intercept (at q = 0) of the two-particle correlation func¬ 
tion, will be considerably different as compared to the 
case of a normal (i.e., in the limit of Q —> 1) pion, but 
will be recovered in the appropriate limit. 

For deriving the relations that allow to describe the 
single- and two-particle inclusive distributions, we start 
by assuming that the Q-boson creation operator in coor¬ 
dinate space can be expressed by [6] 

I 

where a| is the creation operator for creating a Q-boson 
in a quantum state characterized by a quantum number 
1. Then, is one of eigenfunctions belonging to a 

localized complete set, which satisfies the orthonormality 
condition 

J = Sij^, ( 2 ) 

and completeness relation 

='^(x-y). (3) 

I 

Similarly, the Q-boson annihilation operator in coor¬ 
dinate space can be written as 

I 

In momentum space, the corresponding Q-boson cre¬ 
ation operator, V>i(p), and annihilation operator, ^^(p), 
can be expressed, respectively, as 

^^(P) = 

I 

and 

^(p) = 

i 

where 

^i(p) = J V’i(x)e*P'’"dx. (7) 

The Q-bosons are then defined by means of the algebra 
satisfied by their creation and annihilation operators, i.e., 
[17] 

aiaj, - = Sij' 

[ai,ai'] = [al,al] = 0 , 

[Ni,aii] = 

[7Vi,aj',] = 

(8) 


2 



Here Q is a (C-number) parameter, assumed to be within 
the interval [— 1 , 1 ], and Ni is the number operator, which 
can be expressed as 






( 1 - 0 *) 


(9) 


It can be easily verified that, for Q = 1, the normal 
bosonic limit is recovered, i.e., the particles then obey 
the regular bosonic commutation relations, as it would 
be expected. 

We write the density matrix operator for our Q-bosonic 
system as 


p = exp 


pi = exp 


-^{H-pN) 


--(Hi-pNi) 


= l[pi, 

I 


( 10 ) 


where 


Hi=EiNi- N = Y,Ni , ( 11 ) 

i i 

are the Hamiltonian and number operators, respectively; 
T is the temperature. 

The corresponding normalization is explicitly included 
in the definition of the expectation value of observables 
as, for instance, for an operator A 


{A) = 


tr{pA} 


tr{p} 

With the above definitions, it is easy to verify that 

tr{Pi) = ^i{n\pi\n)i = - - . \ ^ 

^ l-exp[-A{Ei-p)] 

where 


( 12 ) 


(13) 


In), = ^(al)"|0) ; [n] = ^ 


\/R! 


1-Q 


(14) 


From the above equations, we can compute the expec¬ 
tation values 


(aja,) = 


1 


exp[y(S, - p)] - Q ' 


and 


(ajalaiai) = 


1 -f Q 


(15) 


(16) 


Ai) _ Q 2 | 

Then, the single Q-boson distribution can be written 


as 


1^1 (P) = (^’^(P)^(P)) 

= EE ^r(p)#(p)(a! 


ai> 


(17) 


The expectation value (ajaii) is related to the occu¬ 
pation probability of a single-particle state I, Ni, by the 
following relation 


{ajai') = Spi'Ni . 


(18) 


For a Q'bosonic system in equilibrium at a temperature 
T and chemical potential p, Ni is represented by the 
modified Bose-Einstein distribution 


Ni = 


1 


exp [^{Ei - p)] -Q 


(19) 


By inserting Eq. (18) and (19) into (17), we obtain the 
single-particle spectrum for one Q-boson species as 


^i(p) = E^'^*(p)^'(p) • 


( 20 ) 


Similarly, the two-Q-boson distribution function can 
be written as 

f 2 (Pl,P 2 ) = {'>PHpi)AHp2)iP{Pi)A{P2)) 

= E ^ri(pi)^r2(p2)^i3(pi)^i4(p2) 

= E ^ri(pi)^r2(p2)^i3(pi)^i4(p2) 

{aj^aiAialaiAh^h 
+ {al^^h){alaiAh^h 

-|- (Sjl 0^2 0,3 d,4 ),1 —12=1^=14 

= Piipi)PiiP2) + |^fv,A^;(pi)V(,(P 2 )p 

I 

+ '^'>P!{pi)Ai{p2)Mpi)Mp2) X 

I 

[{alajaiai) -2{alai)A ( 21 ) 

Using Eq.(15) and (16), we finally have 

f2(Pl,P2) = Pl(pi)Pl(p2) + |^fV,A^;(pi)V(,(p2)P 

I 

-'^'>P!{pi)Ai{p2)Mpi)Mp2) X 


(1 - Q) • Ni 


2 exp[A{Ei - p)]+Q 


exp[i(E, - p)] - (52 


( 22 ) 


The two-particle correlation can then be written as 

^ , N -f 2 (Pl,P 2 ) 

C^ 2 (Pi,P 2 ) — - 5 - 7 — 7777 —r — 

Hi(Pi)T’i(P2) 

= 1+ |E^H^'(Pi)I^E^'I^'(P 2 )I^| X 


I I' 
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1,1' 


1 — — Q) 


exp{^{Ei - fi)) + Q 1 

ex.p{^{Ei -/i)) -Q2 J 


(23) 


It is interesting to note that, for Q=l, we regain the 
results in Ref [6]. Moreover, for Q=0, we also get identi¬ 
cal results as shown in the Appendix A of that reference, 
corresponding to classical Boltzmann distribution, for ei¬ 
ther single or multi-modes. Nevertheless, the naively ex¬ 
pected classical limit of C' 2 (pi,P 2 ) ^ 1 is recovered for 
single mode only, independently on the values of pi and 


P2 

In analogy to the common practice in bosonic interfer¬ 
ometry, it is natural to introduce the intercept parameter, 
A, by means of the relation A = C(q = 0, K) — 1) i.e., as 
the intercept of the two-particle correlation function from 
which the classical limit is subtracted (this procedure, 
however, has no relation to the historical interpretation 
of A as a parameter signaling either total chaoticity or 
partial coherence of the emitting source). 

From Eq. (23), it is straightforward to show that, when 
q = Pi — P 2 = 0, and, consequently, K = ^(pi + P 2 ) = 
Pi = P2 


C'2(K,K) = 2-^-=- 

[ (1 - g)^^r(K)^r(K)^z(K)^i(K) • iVf 


exp(i(A; - n)) + Q 

exp{^{Ei - n)) - g2 


(24) 


We see that, for Q = 1, we regain the ideal result 
for the bosonic intercept at the zero momentum differ¬ 
ence, q = Pi — P 2 = 0, i.e., C' 2 (K, K) = 2. On the 
other hand, for Q = 0 we again recover the naively ex¬ 
pected Boltzmann result, C' 2 (K,K) = 1 for single modes 
only, and in particular for pi = p 2 = AT, as discussed 
above. For multi-modes, however, it is not recovered, 
since there is some remnant communication among Q- 
bosons approaching the classical limit in this case, also 
verified in the result presented in the appendix mentioned 
above. 

From Eq.(24) the intercept of the correlation function, 
A, can be immediately identified as 


A(K) = C2(K,K)-1 = 1- 


_ (1-Q) _ 


^^;(K)^r(K) 

I 




f exp(A(F;; - m)) + Q \ 
\exp{^{Ei - ^)) - QV 


(25) 


We see that the intercept decreases with decreasing Q, 
being always smaller than unity for 0 < Q < 1. On the 


other hand, it is interesting to point that the definition 
for the intercept parameter given by Eq. (25) differs from 
the one in Ref. [15,17], mainly, but not only, because it is 
there defined exclusively for a single mode. In that refer¬ 
ence, comparison is made with experimental points for A, 
which has always been in the limit 0 < A < 1. We could 
proceed similarly within our model as well, by compar¬ 
ing Eq. (25) to the experimental points. Nevertheless, 
we prefer not to do so because it is well-known that other 
factors, such as resonances, dynamical and multi-particle 
effects, as well as kinematical cuts, could also cause the 
intercept to drop into that interval. 


The above derivation can also be reformulated within 
the Wigner function approach. Eor doing this, we de¬ 
velop the product of four wave-functions in Eq. (23) 
into the product of the corresponding Fourier transforms, 
leading to 


^r(pi)^r'(p2)^/(p2)#(pi) = 




n 


(27r)3/2' 

d^r'. 


-P--^V'r(ri) 


d^ 


r2 _ 


(27r)3/2 


2 p 2 .r 2 




(27r)3/2 


f d^r' 


= /A.--- / ^ 


X / d^y y 


d^Ay 

(^ 


AK.Ay / ♦ 


Ay 


Ay. 


V’My-—)#(R-+—)> 

(26) 


where we have defined K = (pi -f P2)/2 as the two-Q- 
boson average momentum, and q = pi — P 2 as their rel¬ 
ative momentum. For writing the last equality, we have 
also changed variables as follows: ri — r 2 = Ax ; ri -|- 
r 2 = 2x ; r'l - r '2 = Ay ; r'^ -H r '2 = 2y. 

Then we can define the Wigner function associated to 
the state I as 


5 /(x,K) 


d^Ax 


Ay 



(27) 


We can proceed analogously to define the equivalent 
function for the integration in y and Ay, remembering 
that ( 7 ;(x, K) = ( 7 ;*(x, K). Then, denoting by 

5(x,K) 91, 

i 


we can finally define the generalized Wigner function of 
the problem as 


ff(x, K; y, K) = ^(x, K)g{y, K) 
exp(A(E; - g)) + Q 


E 


exp(y(Ai - y)) - g2j 


- ( 1 -Q) 

5 /(x, K) 5 ;(y,K)| . 
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(28) 



We see that, for Q = 1) the above expression is reduced 
to the usual result of the original Wigner function, i.e., 
(/(x, K;y,K) = g(x, K)( 7 (y, K). On the other hand, for 
(5 = 0, Eq. (28) is identically zero for single modes only, 
as it would be expect in the limit of Boltzmann statistics. 
Nevertheless, in the multi-mode case, as already shown in 
Eq. (24), there seems to be some sort of residual correla¬ 
tion among (5-bosons even in the classical limit. Aiming 
at better understanding this limit we further explore the 
(5^0 case in the Appendix I. 

By means of this Wigner function, the two-(5-boson 
correlation function can be rewritten as 


C'2(Pl,P2) = 1 + 


//e--q(x-y)g(x,K;y,K)dxdy 
f g(x,pi)dxf g(y,p 2 )dy 


(29) 


The above generalized Wigner function, g(x,K), can 
be interpreted as the probability of hnding a (5-boson 
at a point x with momentum K. Differently from previ¬ 
ous formulations [18-20] we see that, if pions are treated 
as (5-bosons under certain regimes, there is now an ad¬ 
ditional term in Eq. (28). The modified two-particle 
Wigner function no longer can be reduced to the Fourier 
transform of the product of two single-particle Wigner 
functions, but acquires an extra term depending on Q 
in a non-trivial way. As a consequence, for 0 < (5 < 1, 
we can anticipate that the correlation function would be 
narrower and the intercept parameter. A, would drop be¬ 
low unity. We will illustrate more clearly the effects of 
the deformation parameter, Q, on the correlation func¬ 
tion and on the intercept parameter in the next section, 
by means of two toy models. 

In summary, we could say that, maybe under cer¬ 
tain circumstances, pions produced in heavy-ion colli¬ 
sions could be treated as free particles. Nevertheless, and 
as motivated in the beginning of the present section, in 
many others, the interactions of pions among themselves 
and with other particles produced in relativistic heavy- 
ion collisions may not be negligible. In these cases, simi¬ 
lar to what has been suggested in Ref. [17,15], what it is 
proposed here is to mimic those interactions by consider¬ 
ing pions as (5-bosons. In particular, the interpretation 
of Q as an effective parameter reflecting the fermionic 
constituents [16] of the (5-bosons is appealing. Mainly if 
we consider that unconfined degrees of freedom could be 
produced in high energy heavy ion collisions and mani¬ 
fest themselves as (5-bosons in the pre-bosonic stages In 
this sense, they would be regarded as memory traces from 
those pre-confined stages, just before the boson emission. 


III. TWO-(5-BOSON CORRELATION FROM A 
FINITE VOLUME 


A. Toy model 


To explore the effects of the deformation parameter 
Q and of the boundary on the single- and two-(5-boson 
distribution functions, we assume that they are conhned 
in a one-dimension box, [0, L], for simplicity. The three- 
dimensional extension should be straightforward. It can 
be easily checked that the corresponding wave-function 
in the 1-D case is given by 


with 


ipk{x) 



X, 


(30) 


k • L = mr, n = 1,2,3... 


(31) 


Then, the corresponding Fourier transform, tpk{p)t can 
be expressed as 


i’kip) = 


1 1 exp{i{k — p)L) — 1 


(27r)i/2 p-k 

exp{-i{k +p)L) - 1^ 


(32) 


or, equivalently, its square modulus would be written as 

1 ■ 2 (p-k)L ■ 2 (p+k)L 

I r / m 2 1 fSin n ^ Sin ^ 

\*M\ =—[ ,, + (p + tp 

-_P_ _ 

- —^cos(fcL)]. (33) 


ttL {p — kY 
2 sin (—2 


-b 


fcp)A)2sin(^^=^ 


(p-k){p + k) 


On the other hand, if we recall the dehnition of the 
delta function 


,,, , r lsin(a;-L) 
o[x) = lim-, 

L^oo TT X 

it is easily verified that, when L —> oo, we have 

1 


PYp) = ^N,= ^ 


27r exp{\^)-Q' 


(34) 


(35) 


That is, in the limit of an inhnite 1-D box, we obtain a 
modified Bose-Einstein distribution, where the deforma¬ 
tion parameter Q replaces the unity factor, characteristic 
of BE statistics. In the Hnite box case, however, the spec¬ 
trum should change more drastically, due to quantum ef¬ 
fects, which we already showed in Ref. [6]. Moreover, 
in the present case, we have both the finite size and the 
deformation parameter effects combined. To illustrate 
this, we show in Fig.l, the single spectrum distribution 
for two different box sizes. In that plot, as in all the 
others that follow, we have chosen a null chemical po¬ 
tential, i.e., p = 0, for simplicity. For comparison, the 
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Bose-Einstein spectrum distribution, as well as the cor¬ 
responding modified form given by Eq. (35), are also 
shown. It is interesting to note that for finite systems 
and decreasing values of the Q parameter, the width of 
single Q-boson distribution becomes broader, causing the 
maximum of the distribution to drop and the tail to rise, 
due to the conservation of the number of particles. The 
drop of the maximum for the same value of the momen¬ 
tum but for a smaller value of Q would correspond to a 
weaker bosonic behavior of the particles when compared 
to the Q —> 1 limit, leading to a lower occupancy for 
small values of the momenta. The effect is more pro¬ 
nounced for increasing size of the emission region. On 
the other hand, decreasing the values of the deforma¬ 
tion parameter Q has a similar effect as to decreasing 
the source emission size (see Ref. [6]), which is consis¬ 
tent with the uncertainty principle since, as the volume 
of the system decreases, the uncertainty in the pion coor¬ 
dinate decreases accordingly, causing larger fluctuations 
in the pion momentum distribution, which then becomes 
broader. 

It is interesting to check how our result would com¬ 
pare with the interpretation given in Ref. [16], for which 
Q could be viewed as an effective parameter reflecting 
the internal degrees of freedom of the bosons. In that 
reference, the deformation parameter Q is related to the 
ratio of the bosonic volume (L^) to the system volume 
(V) by « 1 — L^/V, where L is the boson’s RMS ra¬ 
dius. The ratio is then correlated to the degree of bosonic 
overlap. Although we do not consider here the bosons 
as extended objects, we still could try and see if that 
picture is compatible with our study. Let us first con¬ 
sider that the bosons have a fixed size. We then compare 
the above relation for two values of the system volume 
(where V 2 > Ri), associating a value of Q to each case. 

It is very simple to see that Q 2 = ^(1 ~ Qi)^ i-®-) 

an increase in the volume would result in a smaller de¬ 
formation parameter, reflecting a smaller overlapping of 
the bosons and their constituents. In other words, for 
a fixed bosonic size and if we enlarge the volume that 
contains the bosons, the resolution decreases, implying 
that Q increases, i.e., gets closer to the boson statistics 
case for which Q = 1. Let us take another approach, by 
considering Q fixed and studying what happens for in¬ 
creasing volumes. In this case, a system of Q-bosons in 
a volume Vi would be associated to a Lf/Vi and another 
one, in similar conditions but with a volume V 2 > Ri, 
would have L 2 /R 2 . In order to keep Q the same, the ra¬ 
tio has to be kept the same, which means that L 2 > Li. 
This could be interpreted as if we had higher resolution 
of the internal degrees of freedom in the second case (i.e., 
the boson with bigger L would have the effect of its in¬ 
ternal constituents more sharply probed). Consequently, 
for the same value of Q, we would expect that the larger 
the system is, more sensitive it would be to fixed value 


of Q. This is precisely what we can see in Eig.l, since 
the effect is more pronounced for L = 8 fm than it is for 
L = 4 fm. 

We have seen that the deformation parameter has a 
significant effect on the spectrum of the bosons. We dis¬ 
cuss next what this implies to the interferometry of two- 
identical Q-bosons. In Fig. 2, the correlation functions 
for two values of the mean momentum K are shown for 
different deformed parameter, Q, as a function of the 
pair relative momentum, jqj. For Q = 1, as already 
shown in Ref. [6], we see that, as the mean momen¬ 
tum increases, the source radius increases accordingly, 
due to the fact that contributions from small momenta 
come from smaller quantum states I which, in turn, cor¬ 
responds to larger spread in coordinate space. Similar 
behavior in the radius is seen as Q decreases below unity. 

Another interesting point concerns the way the param¬ 
eter Q changes the intercept parameter. A, of the two-Q- 
boson correlation function. The effects on A are more 
pronounced as smaller values of K are considered, which 
is natural, as for large values of the average momentum, 
the quantum effects become less relevant. From Eq.(25), 
we can see that, for increasing K, the dominant factors 
come from the larger I states which, on the other hand, 
give smaller contribution to the two-Q-boson correlation, 
due to the factor N\, which decreases for increasing AT, 
as can be seen from Eq. (9). Consequently, this makes 
the intercept parameter to vary more slowly for increas¬ 
ing values of K. This is illustrated in Fig. 3, where A 
is shown as a function of Q for different values of the 
mean momentum and for different source radii. Again, 
we note that as the source radius becomes bigger, the 
Q effects on the intercept parameter become less signif¬ 
icant, since in this case the quantum effects are smaller. 
In the plot, we only shown the variation of A for Q in the 
interval [0,1], corresponding to A < 1 . Of course, if Q 
is larger than one, as one could expected from Eq.(25), 
the value of A could be bigger than one. Also if the value 
of Q is negative, A could be less than zero. However, we 
are treating here bosons with a modified commutation 
relation. Since such unexpected behavior for the inter¬ 
cept parameter was never observed experimentally, for 
any type of bosons, we do not consider this case here. In 
other words, our analysis refers basically to the interval, 
0 < <5 < 1. Nevertheless, we should keep in mind that 
the intercept parameter of the two identically charged 
Q-bosons could be bigger than one or less than zero, for 
some specific values of Q. 

Although not shown in Fig. 3, the limit Q ^ 0 de¬ 
serves a closer analysis. As we briefly discussed in Section 
II, the naive classical limit C' 2 (pi,P 2 ) ^ 1 is recovered 
only for single modes. The multi-mode case is analyzed in 
detail in Appendix I but we summarize the main results 
here. We start with the limit Q ^ 0, for which the first 
of the commutation relations for equal modes in Eq. (8) is 
reduced to aa+ = 1. Nevertheless, we demonstrate in the 
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Appendix I that it is not only the commutation relations 
that matter when discussing the behavior of the inter¬ 
cept parameter, A. The density matrix seems to play an 
essential role, at least for the type of Q-boson we analyze 
here. In this case, with the density as defined in Eq.(ll), 
we get A ^ 0 for very small system sizes (and Q = 0), 
recovering what would be naively expected for classical 
particles. In the opposite limit, i.e., for very large sys¬ 
tems, we get a constant Nk, resulting in A ^ 1/3. 

The above limits could suggest that we get different 
results depending on the wave function and boundary 
conditions, reflecting the dependence on the dimensions 
of the system and, consequently, on Ek, since the par¬ 
ticular pk we chose contains an explicit dependence on 
the energy of the state. However, we also demonstrate 
in Appendix I that, if we had chosen a different density 
matrix than the one in Eq.(ll), for instance, pk = const., 
the wave function would not play any role, since we have 
Nk =const., and there is no energy dependence in this 
case. As a result, we get A ^ 1/3, independently on the 
system size and Ek- Although we do not show in Fig. 3 
the limit A (AT = 0), we can still check the consistency of 
the results plotted there with the analysis we have just 
made. For that, we should look into the smallest value 
of momentum shown in that plot, i.e., K = 0.3 GeV/c. 
Then, according to our analysis in the Appendix I, the 
value of the intercept parameter for small systems would 
tend to approach the limit A ^ 0, whereas for large ones, 
it should approach A ^ 1/3. Adapted to our plot, this 
would mean that \{L = 4/m) < X{L = 8fm) for Q = 0 
and K = 0.3 GeV/c, which is precisely what is seen in 
Fig. 3. We can also verify that this result is more gen¬ 
eral, since the same feature is again reproduced in Fig. 
6, as we shall see in the next subsection. 

Furthermore, if we return to the expression of the gen¬ 
eralized Wigner function, Eq. (28), we see that the 
usual Wigner function is recovered (last term vanishes) 
for (5 = 0 but, again, for single modes only. For multi¬ 
modes, however, the Wigner function is modified, even 
for (5 = 0. This result seems to indicate the existence of 
some kind of residual correlations among the particles in 
the system i.e., an extra communication among the par¬ 
ticles of different states, beside the commutation relation 
defined among particles in the same state. 

Along the lines discussed above, only recently we be¬ 
came aware of a Monte Garlo event generator by Wilk et 
ah, which makes an attempt to improve Bose-Einstein 
correlations in numerical modeling. It would be very 
interesting to run their simulation and check for con¬ 
sistency with our numerical calculation of the correla¬ 
tion function, C 2 ipi,P 2 ) and of the intercept parameter, 
X{K), mainly regarding the multi-particle effects. 


B. (5-bosons are confined inside a sphere 

In this section, we consider the that the pions pro¬ 
duced in high energy heavy-ion collisions, treated here as 
the hypothetical (5-bosons, could be bounded in a sphere 
up to the time immediately preceding the freeze-out of 
the system. This is conceived in such a way that their 
distribution functions are essentially the ones they had 
while confined. Analogously to the procedure developed 
in [6], the pion wave function in this case should be de¬ 
termined by the solution of the Klein-Gordon equation, 
i.e.. 


[A -I- tp{r) = 0 , (36) 

where = E'^ — m? is the momentum of the pion. On 
writing the above equation, we have assumed confine¬ 
ment, i.e., the potential felt by the pion inside the sphere 
is zero, while outside it is infinite. The boundary condi¬ 
tion to be respected by the solution is 

^{v)\r=R=Q , (37) 

where R is the radius of the sphere at freeze-out time. 

The normalized wave function corresponding to the so¬ 
lution of the above equation can easily be written as 

= 0 (r > R). 

(38) 

The momentum of the bounded particle, k, can be de¬ 
termined as the solution of the equation 


Ji+i(/cA) = 0. (39) 

Inserting Eq. (38) into Eq. (7), we can determine the 
Fourier transform of the confined solution of a pion inside 
the sphere, as a function of the momentum p, as [6] 


'^fc/m(p) — a/ i X^lmip) 

P 


-k 


p2 _ J^2 


Ji+UpR) . (40) 


In terms of Eq. (19) and (20), the single-inclusive dis¬ 
tribution function can be written as 


^I(P) 


kin 


= E 


2;-f lA f kJi+NpR)' 


-/ exp -Q \ 2Trp J y p'^ - k^ 


(41) 


In the limit R oo, the single particle spectrum, can 
then written as 
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1 


V 


^l(p) 


exp - Q 


(27r)3 


(42) 


where V = is the volume of the sphere. We see from 
Eq. (42) that the modified Bose-Einstein distribution is 
recovered in the limit of very large volumes. 

In Fig. 4, the normalized single-particle distribution 
is plotted as a function of |p|. We clearly see that, due 
to the boundary effects, the maximum value of |p| in the 
spectrum decreases for decreasing volumes, being always 
smaller than the case corresponding to the R oo limit. 
On the other hand, this is similar to the result obtained 
with the previous example of the 1-D box, and the con¬ 
finement does not seem to cause a significant effect on 
the spectrum. 

We can write the expectation value of the product of 
two Q-boson creation operators in momentum space as 
before, resulting in 


w'(p.)dP=)> = E^M|(|Y^ 


klm 6Xp 


V ^ 


Q 


= E' 


klm exp 


fMj^] -Q 



■ -k 

_p\ — k'^ 

V Pi 


FT, , 

\ 

-k 1 

V P2 

_pI - 


Ji+i{piR) 

Jl+liP2R) 


= E' 


1 


ki exp 


( Mmz 


V T 


-Q 


PlP2 


(PiR) 


Jl+i{p2R) PliPl ■ P2) ■ 


(43) 


The two-pion interferometry correlation function can 
then be estimated by inserting the expressions (41) and 
(43), into Eq. (23). We see from the above results that, 
in general, this function depends on the angle between pi 
and p 2 , similarly to what was discussed in [6]. For the 
sake of simplicity, however, we will also consider here pi 
parallel to P 2 , implying that Pi{pi ■ P 2 = ±1) = (±1)*- 
The results for two-pion interferometry corresponding 
to two different values of the pair average momentum 
K = (pi -l-p 2 )/ 2 , but fixed temperature, are shown in 
Fig. 5. For Q = 1 we can see that, as the pair aver¬ 
age momentum, K, increases, the apparent source radius 
becomes bigger, which reproduces the result obtained in 
Ref. [6]. However, considering K = 0.3 GeV/c, if we 
compare the cases corresponding to Q = 1 and Q = 0.5, 
respectively, we see that the resulting correlation func¬ 
tion becomes narrower and the intercept parameter. A, 
drops below its previous unit value. On the other hand. 


if we now keep this value of Q = 0-5 but consider K = 0.5 
GeV/c, the width is even narrower, but the intercept is 
higher than that corresponding to K = 0.3 GeV/c, al¬ 
though it still is below one. This due to the fact that 
for smaller momentum pairs, the quantum effects are 
stronger, as in the previous toy model studied in Fig. 
2. The effect comes from the contribution of the third 
term in two-Q-boson interferometry formula, in Eq. (23). 
We could understand these results by noting that pi- 
ons with larger momentum come from larger quantum 
I states which, in turn, correspond to a smaller spread in 
coordinate space. Due to the weight factor in Eq. (23), 
of modified Bose-Einstein form, larger quantum states 
will give a smaller contribution to the source distribu¬ 
tion, causing the effective source radius to appear larger. 
On the other hand, this behavior is interesting if we com¬ 
pare to results corresponding to expanding systems. In 
this last case, the probed part of the system decreases 
with increasing average momentum [22,23]. Naturally, 
our present approach does not consider the effects of ex¬ 
pansion and the enlargement of the system’s apparent di¬ 
mensions with increasing K, seen in Figure 5, has its ori¬ 
gin in the strong sensitivity to the dynamical matrix. In 
Ref. [24], the combined effects of a finite boundary and an 
expanding system were considered together. What they 
observed was an opposite effect, i.e., the effective source 
radius extracted from two-pion interferometry would de¬ 
crease as K increases. However, for small values of the 
momentum K, it seems that the boundary effects are 
dominant over the expansion effects. 

In Fig. 6, we plot the intercept parameter A vs. Q 
for different values of momenta and source radii. We see 
that the similarity to the results in Fig. 3 is evident: 
A becomes bigger as either the total momentum of the 
pair, K, or the source radii, R, increases, reaching val¬ 
ues significantly below one for sufficiently small values 
of either one of those variables. This is expected, since 
the quantum effects are more prominent for smaller mo¬ 
menta. Again, the comments relating these results to 
the speculation in Ref. [16] apply here, as in the previous 
example. 

In summary, by comparing the interferometric results 
corresponding to the two toy models, we see that a gen¬ 
eral behavior is roughly reproduced in both cases. First, 
Figures 1 and 4 show that the results for the normal¬ 
ized spectra are very similar, the maximum of the curves 
dropping with decreasing values of Q, followed by the rise 
of the respective tails, this being related to the conserva¬ 
tion of the number of particles. Second, the width of the 
correlation functions, as seen in Figures 2 and 5, seems to 
decrease for increasing values of Q. Besides, we see from 
the Figures 3 and 6 that the Q-boson effects cause the 
intercept parameter. A, to drop below unity. This clearly 
demonstrates that treating pions as Q-boson indeed alter 
the two-pion interferometry. Besides, we saw in our ex¬ 
amples that it also modifies and generalizes the form of 



two particles Wigner function. Nevertheless, we also saw 
that the effects for the bounded case are less pronounced 
as compared to the unbounded one, at least for the set 
of parameters adopted in the calculation. 

IV. CONCLUSIONS 

In this paper, we derive spectrum and correlation func¬ 
tion relations by adopting the density matrix given in 
Eq.(lO), suitable for describing charged Q-bosons. The 
finite volume effects on the Q-boson spectrum were then 
studied in Figures 1 and 4, for two specific examples, 
leading to similar results as in Ref. [6,8,10,11], for Q = 1. 
We find that the small momentum region is depleted for 
the modified Bose-Einstein distribution with respect to 
case when Q ^ 1. The effects on two equally charged Q- 
boson correlation function were also analyzed here. The 
results in Figures 2 and 5 show that the correlation func¬ 
tion shrinks for increasing average pair momentum, cor¬ 
responding to an increase of its inverse width [6]. We also 
observe that its intercept drops for decreasing Q. In other 
words, it was shown in Figures 3 and 6 that the intercept 
parameter. A, decreases for increasing Q, for the same 
values of the average transverse momentum and source 
radius. On the other hand, A becomes larger when either 
the mean momentum of two-Q-boson or the source radius 
increases. This result reflects a strong sensitivity to the 
dynamical matrix, through the modified Bose-Einstein 
weight factor. For Q = 1, previous pion interferometry 
results are regained, as well as the Boltzmann-type dis¬ 
tribution for Q = 0, as can be verified in Ref. [6]. 

We have also derived a generalized version of the two- 
boson Wigner function, allowing to treat the case of two 
Q-bosons. This result is particularly important because 
it shows that, in more general situations, the Wigner 
function is distorted with respect to its decoupled form, 
i.e., 5 (x,K;y,K) = ^(x,K)g(y,K). The extra terms it 
acquires, as shown in Eq. (28) (see also Eq. (70) in the 
Appendix II), reflect non-trivial interactions among the 
Q-bosons. This generalization, however, is reduced to its 
well-known form in the limit Q ^ 1, as it should. 

We also analyzed in the current work how our results 
compare with the interpretation given in Ref. [16], for 
which Q could be viewed as an effective parameter re¬ 
flecting the internal degrees of freedom of the bosons. If 
we consider a fixed value for the deformation parame¬ 
ter, Q, and compare results for increasing volumes, our 
result would be compatible with that interpretation. To 
summarize the comparison we could say that, for increas¬ 
ing volumes and keeping the same value of Q, we would 
have to increase the resolution, i.e., the sensitivity of the 
probe to the internal degrees of freedom of the boson. 
Consequently, for the same value of Qj we would expect 
that the larger the system is, more sensitive it would be 
to this parameter. This is precisely what we can see in 


Fig.l, since the effect is more pronounced for L = 8fm 
than it is for L = 4fm. 

The derivation analyzed here has several common fea¬ 
tures with the one in Ref. [15,17] but here we adopt a 
entirely different approach, focusing in what seems to us 
the most important part of an interferometric analysis, 
i.e., the correlation functions themselves. We also ana¬ 
lyze the effects on the spectra and on the intercept pa¬ 
rameter, A, but, again, our result is more general than in 
that reference, since it is there restricted to single modes. 

Another remark concerns the relation of the parameter 
A and its possible interpretation as a partial coherence 
of the emitting source for values below unity, as com¬ 
monly found in the literature on boson interferometry. 
It is well-known that many effects, such as resonances, 
dynamical and multi-particle effects, as well as kinemati- 
cal cuts, could also contribute to yield values of so-called 
chaoticity or coherence parameter below unity, although 
those effects have no relation to partial coherence of the 
source. Similarly, the behavior we discussed of the inter¬ 
cept parameter A, when considering decreasing values of 
the deformation Q, are not related to the way the source 
emit those bosons. The effects discussed here are meant 
to show that they also could cause a deviation from the 
idealized picture. Moreover, for keeping the analysis sim¬ 
ple, we are not taking all the above effects into account 
in our study. For this reason, we preferred to not com¬ 
pare our results with experimental data, and thus did not 
introduce any quantitative analysis of this picture. 

Also, for completeness, we derive in the Appendix II, 
the equations associated to the so-called type-B Q-boson 
interferometry [15,17], corresponding to slightly differ¬ 
ent commutation relations. The derived relations for the 
single- and two-particle distributions are in Eq.(68), and 
Eq. (69), respectively. Also for this case, we propose a 
generalized form for the Wigner function, as can be seen 
in Eq.(70). 
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p(GeV/c) 


FIG. 1. (Color online) The normalized spectrum (in arbi¬ 
trary units) vs. momentum |p| (in GeV/c) is shown. The 
input temperature is T = 0.14 GeV and the chemical po¬ 
tential is /r = 0. The solid line corresponds to the modified 
Bose-Einstein distribution, i.e., to the limit R —> oo. The 
dotted and dashed lines correspond, respectively, to L = 8 fm 
and L = 4 fm. The thicker lines refer to Q = 1.0 and the 
thinner ones to Q = 0.5. 



q(GeV/c) 


FIG. 2. (Color online) Two-pion correlation vs. momen¬ 
tum difference |q| (in GeV/c). The input temperature is 
T = 0.14 GeV and the chemical potential is /r = 0. The 
solid line corresponds to mean momentum K = 0.3 GeV/c 
and dashed one to IF = 0.5 GeV/c. The thicker lines refer to 
the case of Q = 1.0 and the thinner one to Q = 0.5. The box 
size is L = 4 fm. 



FIG. 3. (Color online) The intercept parameter. A, is shown 
vs. Q, the deformation parameter. The input temperature is 
T = 0.14 GeV and the chemical potential is p = 0. The solid 
line corresponds to mean momentum K = 0.3 GeV/c and 
dashed one to the case K — 0.5 GeV/c. The thicker lines 
refer to the case L = 4 fm and the thinner ones to L = 8 fm. 



p(GeV/c) 


FIG. 4. (Color online) The normalized spectrum (in arbi¬ 
trary units) vs. momentum |p| (in GeV/c) is shown. The 
input temperature is T = 0.12 GeV and the chemical poten¬ 
tial is = 0. The solid line corresponds to the case i? = 3 
fm, the dotted one to the case R — 6 fm, and the dashed line 
corresponds to the case, i? = oo. The thicker lines refer to 
Q = 1.0 and the thinner ones toQ = 0.5. 
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FIG. 5. (Color online) Two-pion correlation vs. momen¬ 
tum difference |q| (in GeV/c) is shown. The input tempera¬ 
ture is r = 0.12 GeV and the chemical potential is /r = 0. The 
solid line corresponds to mean momentum K = 0.3 GeV/c 
and the dashed one, to the case K = 0.5 GeV/c. The thicker 
lines refer to the case of Q = 1.0 and the thinner ones to 
Q = 0.5. The sphere size is 7? = 3 fm. 


APPENDIX I 


We here analyze the Q — Q limit in more detail. We 
can see its implications on the expectation values, di¬ 
rectly from Eq.(15)-(16), on the correlation function, 
from Eq.(23), and on the intercept parameter. A, from 
Eq.(25), by simply imposing the limit (5 = 0. Alterna¬ 
tively, we can start with our definitions in Eq.(8), and 
thus simultaneously test our results. Lets also go back to 
the definition in Eq.(14). For Q=0 we see that 

Thus, [n]! = 1 in both the above cases. Consequently, 
instead of Eq.(14), we have 

In), = (al)"|0) . (45) 

If we apply to the above equation the annihilation and 
creation operators, we have 

ai\n), = ai(a;^)”|0) = |n - 1) ; al\n)^ = |n -h 1) , (46) 



FIG. 6. (Color online) The intercept parameter. A, is shown 
vs. the deformation Q. The input temperature is T = 0.12 
GeV and the chemical potential is p = 0. The solid line 
corresponds to mean momentum K = 0.3 GeV/c and dashed 
one to the case K = 0.5 GeV/c. The thicker lines refer to the 
case 77 = 3 fm and the thinner ones to 7? = 6 fm. 


which follow from the first commutation relations in Eq. 
(8) and from the condition a;|0) = 0. The other commu¬ 
tation relations remain unchanged in this limit. 

Let us now estimate the trace of the density matrix, 
tr[pi). From Eq.(lO) and (11), we have 


CO OO 

n n—0 m—0 


1 

m! 


OO 

n—0 


1 

1 - exTp[-^{Ei - p)] 


(47) 


For estimating the limits of Eq.(15) and (16) for Q ^ 0 
along the lines above, we have to estimate the traces 
tr{a\aipi) and tr{a\a\aiaipi), as follows 


tr{a]aipi) = '^^{n\a\aipi\n)^ = 

n 

/it r ~ m) im \ 

n—0 m—0 

= E E[-^(^' - - l|n - 1), 

I ml 




n—0 




e T 


1 - exp[-^(E;; - p,)] 


(48) 
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tr{ajalaiaipi) = ^ ,(n|aj'a|'a/a/pi|n), = 

n 

1 AT 

n—0 m—0 


LXJ (_XJ ^ ^ 

= XI XI “ l|ajai|n - 1), = 

n—1 m—0 

OO OO 1 1 

“Si; 1-r*-®' - »‘>"l’";n.<" - 2I" - 2), = 

n—2 m—0 


= ^ g- + (-Ei-/i)n _ g-i(Ei-M) _ X = 


n=0 


g-TC^Z-Ai)] 

— g~ (-^z ~Az) 


(49) 


From the above results and our definition in Eq.(12), 
it follows immediately that 


ajai) = + = Ni , 

(50) 

{ajajaiai) = _ 

(51) 


It is then straightforward to conclude that both results 
coincide with the ones in Eq.(15) and (16), reinforcing the 
correctness of our results. 

The purpose of the above derivation goes beyond what 
we have just discussed, it is also helpful for discussing 
more deeply the peculiar limit of Q —> 0. As we have 
already pointed out in the body of the manuscript, we 
recover the expected classical result in that limit for sin¬ 
gle modes only. In other words, only for single modes 
we have, for g = pi — P 2 = 0 (and, consequently, 
K =pi = P 2 ) 

lim C 2 (K, K) = 2 ; lim XlR) = 0. 

Q^O ^ ' Q^O ^ 

In this case, we also recover the well-known form of the 
Wigner function, i.e. 


from which we see that only odd values of n, i.e., n = 
2m — 1 can contribute, resulting in 


l^fc(0)l 



( 2 m-l) 2 ^^ (m=l,2,...) 
0 (n even) 


(53) 


Then, for writing the chaoticity parameter in the limit 
of Q = 0 for the toy model, we rewrite Eq.(25) for this 
particular case, as 


A(K) = 1 - 




(54) 


We have then to estimate the sums separately, as fol¬ 
lows 

XiV/|^z(K)p=^e-^|^z(K)|E (55) 

i i 


X \MK)\^ = X l^/(K)|" . (56) 

i i 

For estimating the above sums, it is more convenient 
to consider appropriate limits for the size of the 1- 
dimensional box. Let us first assume the limit of very 
small sizes, i.e., L ^ 0. In this case, Ei « mr/L, we 
denote as 

— ^1 ^ . N 

X = e 

A 2 m-i = (2 to — l)7r/L = {2m — l)i?i 

ZV2r„-l = (57) 



Consequently, the sums for K = 0 can be written as 


X^™i^„.(o)p = X^“ 


4L 1 
TT^ {2m — 



m 


^2m-l 

(2m- 1)2 ■ 



(58) 


^im 5 /(x,K) = g{x,K)g{y,K). 

Nevertheless, when multi-modes are taken into ac¬ 
count, the behavior of the above quantities change con¬ 
siderably. For better illustrating this fact, it is convenient 
to adopt a specific example. Let us choose our toy model 
in section III. A, for simplicity, restricting ourselves to the 
Q = 0 case. In general, A = \{q = pi — P 2 = 0) = X{K = 
Pi = P2) but, for simplifying our analytical study, we re¬ 
strict our analysis to A (AT = 0). In this case, we see that 
the square modulus of the wave-function in momentum 
space can be written as 

|2/’fc(0)P = “ cos{kL)]-,kL = mr{n = 1,2...), 

(52) 


X NL-1 \^2m-l{0)\^ = J2 


2(2m-l) 


16L2 


16L^ 


X 


(2m — l)47r® 

2,2(2m-l) 


(2m — 1) 


4 • 


(59) 


If we then bring Eq.(58) and Eq. (59) into Eq. (54), 
we get 


lim A(K = 0) = 1 — 


E 


m (2m —1)*^ 


L—^0(ai—»-0) 


= 1 — lim 


\2—im (2m.—1)2 


x'^{l 


+ 1^ +V + -) 


X I X 

W 5- 


x^l + 2§; + ^ + 2^ + ...) 


= 1-1 = 0 , 

(60) 
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and thus, the expected limit for classical particles is re¬ 
covered for (5 = 0, when L ^ 0. 

In the opposite limit of L ^ oo, however, we can see 
that k = (2to— l)7r/L ^ 0, which implies that Ek TO,r- 
Consequently, Nk ^ Nq = const. In this case 


lim A(K = 0) 

L—»0(a:—»0) 


1 


I 


2-^m 77^ {2m— TF 




Aro4L 


m 7r^(2m—1)2 



1 

3 ■ 


(61) 


From the result in Eq.(60) we see that, for the par¬ 
ticular density matrix we chose in Eq.(ll), the inter¬ 
cept parameter, A(K = 0), can still have the expected 
value for bosonic interferometry in the classical limit, 
i.e., A(K = 0) = 0, even for multi-modes, due to the 
combined effects of the density of states, Nk, and of the 
wave-functions, = 0). However, in the second case, 
we see from Eq.(61), that limi^oo A(K = 0) = 1/3. This 
is reflecting the relation between L and k coming from 
the shape of the wave function and boundary conditions 
in a finite system, as in Eq.(52). Moreover, the difference 
among the two limits reflects the particular choice of the 
density matrix. In order to demonstrate this, we choose 
a much simpler example than in Eq.(ll). For instance, 
let us choose pk = const., and proceed analogously as be¬ 
fore. In this case, instead of Eq.(50) and (51), we would 
have 


(ajai) = 1 = Ni ; {ajajaiai) = 1 , (62) 


As a consequence, we get for the sums in Eq.(55), (56) 
and Eq.(61) 


Eimor = 


4L 


Ei\M0)\^ = 

16N^L‘^77‘ 


16L* 


A(K = 0) = 1 - 


j2 T 2^4 

967r^ 

/ 4NoL77^ ' 


=1-2.1 
3 3 




(63) 


From what we just saw above, we could say that, if 
we choose a different type of density matrix, i.e., pk = 
const., for instance, the wave function does not seem to 
play any role, since we get A'fc=const. Consequently, for 
this much simpler density matrix, we see that we get 
A ^ 1/3, independently on the energy of the states, Ek, 
and on the size of the 1-dimensional box, L. 


APPENDIX II 


For completeness, we will also derive below the so- 
called type-B (5-boson interferometry formulation, as de¬ 
fined in Ref. [17]. For type-B (5-boson, the operators bi 
and bj satisfy the following commutation relations 

bibl-Q^‘-‘'bl,bi=Si,cQ-^‘ 

bibl-Q-^^.nlbi=SpcQ^‘ 

[bi,bi'] = [bib],] =0, 

[Ni,bi>] = -6pvbi 
[Nubl]=5i,i,b\, 

[Ni,Ni,] = 0. (64) 

In the above relations (5 is a parameter, which can be 
assumed within [—1,1]. If we define it, following Ref. 
[17], as (5 = [0 < 8 < 27r) or, equivalently, cos(6l) = 

^{Q + Q~^), then 

(blk) = 

_ _ exp{^{Ei - p)) - 1. __ 

exp(|;(E; — p)) — 2 cos{9) exp{^{Ei — p) + 1 
^ exp[^{Ei - p)] sinh[^(E; - p)] 
cosh[y(if/ — p)] — cos{d) 

and 

{blblkbi) = 

2cos(0) 

exp(|;(if; — p)) — 2cos(20)exp(A(if; — ^) -I-1 
^ cos(6l) exp[^{Ei- p)] 
cosh[^(ifi — p)] — cos(26l) 

As before, the expectation value {blbc) is related to 
the occupation probability of the single-particle state /, 
by a similar relation, i.e., 

{bybi,) = 6pcNl^'> . (67) 

Then similar to the derivation of Eq.(16) and 
Eq.(20),we have 

^i(P) = ^iblbi)^iiP)Mp) (68) 


and 

-f2(Pl,P2) = C’l(Pl)C’l(P2) + I X!^?'(Pi)^'(P2)I^ 

I 

I 
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= Pi{Pi)Pi{P 2)+Y^ NiNi4!l{pi)'tjjp{p2)'tjji{p2)'ipi'{Pi) 
1,1' 


1-V(i 


008 0) 


4cosh^(%^) + 


COS 0(cos 9 — 1 ) 

/ -Ai S 


sinh2(i 


cosh( ^ ) — cos(20) 



Analogously to what was done at the end of Section II, 
we can also define a modified Wigner function for type-B 
Q-boson. The new Wigner function for this case can be 
defined similarly as before, resulting in 


(x, K; y, K) = (x, K)g(*') (y, K) - 

[4cosh2(%ii) + 


(ivf^)2(l-cos6i) 


sinh2(- 


cosh( ^ ) — cos(26l) 


>.(70) 


Then the two-pion interferometry formula follows anal¬ 
ogously to Eq.(29) for this case, i.e., 

^ . , //e-q(-y)gW(x,K;y,K)dxdy 

^ + fg0>)(x,p,)dxfg0>)(y,p2)dy 

(71) 
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